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Abstract 

Using the quark-antiquark interactions obtained in the framework of 
the bootstrap method we construct a potential model, investigate the 
possibility of describing of heavy quarkonia and calculate the bot- 
tomonium spectrum. The potential of the interaction was obtained as 
a nonrelativistic limit of the relativistic quark-antiquark amplitudes 
QQ -» QQ. 

PACS numbers: 14.40.Lb, 14.40.Ev, 12.39.Pn 



1. Introduction 

Quarkonium systems, the bound states of a heavy quark and antiquark, 
have played a particularly important role in considering of strong interac- 
tion dynamics. The discovery of heavy quarkonia: the families of J/ ip- and 
T-mesons promoted the quark model and non-Abelian gauge field theories, 
leading to the prevailing picture of particle physics. 
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QCD features two remarkable properties. First, asymptotic freedom implies 
that at very high energies and momenta, quarks and gluons interact only 
weakly and act as quasifree particles [1, 2]. Second, confinement presumably 
results from the fact that at low energies the force between quarks increases 
with their distance, so that quarks are always tied into hadrons and cannot 
be removed individually. 

Confinement makes it hard to calculate quantities for the bound states within 
QCD as one cannot apply perturbative QCD. By analogy with the positro- 
nium, and given the large masses of the charm and bottom quarks, nonrel- 
ativistic phenomenological potential models have been applied as tools for 
the quarkonium spectroscopy. To accommodate the properties of QCD these 
models, e.g. [3-6], are based on a short range part motivated by perturbative 
QCD and a phenomenological long range part accounting for confinement. 
In QCD with heavy c— and b— quarks, the characteristic scale Aqcd ~ 0.2 
GeV is small as compared to the quark masses, m c ~ 1.5 GeV and ~ 5 
GeV. A systematic expansion in the powers of I/tuq is possible [7-9]. The 
bound state problem can be dealt with non-relativistically. Following these 
thoughts one treats the bound state problem by solving a Schroedinger equa- 
tion using an appropriate potential. 

In this work the quark interaction potential is considered. The short range 
part of the potential is obtained as a nonrelativistic limit of the relativistic 
quark amplitudes of the bootstrap quark model [10, 11]. These quark ampli- 
tudes depend not only on a squared momentum transfer t, but on the energy 
variable s also. Therefore the direct transition to nonrelativistic potentials is 
not possible: these amplitudes correspond more to quasipotentials [12]. To 
obtain quark potentials from the quark amplitudes one has to fix the energy 
s = so and then the dependence on momentum tranfer at fixed energy is 
considered to be potential. The energy fixing s implies an introduction of a 
momentum cutoff parameter in the Fourier transformation. As a result, 
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the following expression is obtained for the short range part of the potential 
[13]: 

A<i> 

V B (r) = j / - sinkr ~, ^7T^ dk i W 

m 2 J r 1 - gB(k 2 ) 

o 

where g is a dimensionless coupling constant, which is also a parameter of 
the model, B(k 2 ) is the Chew-Mandelstam function [14] for the gluon state: 




where m is the mass of a heavy quark m b , the coefficients /3\ and j3 2 for 1 
state are: f3\ — -,@ 2 — tt- 
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The qualitative behaviour of the bootstrap potential with r is shown in fig.l. 
In difference with the majority of the quark interaction potentials [5, 6, 15- 
20] the bootstrap potential has the finite value at r = in the consequence 
of energy cutoff introduced in calculations of bootstrap quark amplitudes. 
Besides this, the presence of the momentum cutoff leads to the origin of 
small oscillations of the potential at distances ~ 1 fm. 

The potential of confinement is considered as a linear potential with a slope 
defined by the angle a. This potential is added to the bootstrap potential at 
a distance r . a and r are also the parameters of our potential model. So 
the quarkonia potential has the form something like that shown in fig. 2. 
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2. Results 

The quarkonia potential considered in the previous section is used in the 
time-independent Schroedinger equation to find bound states, while spin- 
spin and spin-orbit interactions (Breit-Fermi interaction) are treated within 
perturbation theory. The resulting mass formula is given by 

M{k" s+ H 3 ) = 2m Q + E kl + 3 j^(\s(S+ 1) - ^ |^(0) | 2 + 

j(j + l)-l (l + l)-S(S + l)/l* 
+a s 



frig \ r 3 

where a s is a running coupling constant. 

We use as input to determine the model parameters the states rj b (lS), 
T(15'), T(2S') and C(1P), the center of gravity for the IP triplet states de- 
fined as 

C(1P) = i(5M( Xb2 )+3M(x6i)+M( Xb o)) ~ 9900 MeV. (4) 
Identifying the C(1P) with the 1 1 P 1 state of the model the resulting param- 
eter set is 
Ap = 3.05 GeV, 
9 = 2-22, 

a = 0.28, (5) 
r = 0.39 fm, 

the values of A& (in the 5-function) and nib are taken from the bootstrap 
method [21]. 

The resulting masses are compared with the experimental data in table 1, 
while in fig.3 the resulting spectrum is displayed together with the exper- 
imental one. The bottomonium S-wave state reduced radial wavefunctions 
calculated within the model are shown in fig A and for the P- and D-wave 
states are shown in fig.5&,b. 
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3. Conclusion 

In this work, we have costructed a potential model for heavy quarkonium 
based on the bootstrap potential. The aim of this work has been to investi- 
gate the possibility of describing of heavy quark- ant iquark systems with the 
help of the derived potential using the Schroedinger equation and to get a 
satisfactory description of the quarkonium spectra with minimal phenomeno- 
logical input. 

We have demonstrated that a satisfactory description of the quarkonium 
spectra is possible within this model with reasonable values of the param- 
eters. In the future one should use it to calculate other properties such as 
quarkonia radii, decay widths and branching ratios. 
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Figure and table captions 

Figure 1: The bootstrap potential as a short range part of the model po- 
tential. 

Figure 2: The model potential. 

Figure 3: Bottomonium spectra of experiment and our model; r)b(2S), r)b(3S), 
hb(lP) and h h {2P) mass measurements [23] are included separately in the 
experimental spectrum. 

Figure 4: Bottomonium S-wave state reduced radial wavefunctions calcu- 
lated within our model. 

Figure 5a: Bottomonium P-wave state reduced radial wavefunctions calcu- 
lated within our model. 

Figure 5b: Bottomonium D-wave state reduced radial wavefunctions cal- 
culated within our model. 

Table 1: 66-state masses from the experiment and our model. 
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V(r), GeV 




fm 
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Table 1. 66-state masses from the experiment and our model. 



State 


Candidate 


Experimental mass [22], MeV 


Theoretical mass, MeV 


K - / [J 


ifr(lS) 


9300 ± 20 ± 20 


9300 


1 3 S 1 


TflS") 


9460.30 ± 0.26 


9460 








9900 


1 3 P 


/ A, UU V / 


9859.44 ± 0.42 ± 0.31 


9869 


1 3 P! 


■Xbi(lP) 


9892.78 ± 0.26 ± 0.31 


9884 


1 3 P 2 




9912.21 ± 0.26 ± 0.31 


9916 


2 1 S' 






9997 


2 3 S 1 




10023.26 ± 0.31 


10023 


- 1 - z 






10156 


1 3 P>! 






10150 


1 3 £ 2 


T(1D) 


10161.1 ±0.6± 1.6 


10154 


1 3 D 3 




10160 


2 1 Pi 






10219 


2 3 Pn 

U 


Yhn(2P) 


10232.5 ±0.4 ±0.5 


10191 


2 3 P! 


Yhi (2P) 


10255.46 ±0.22 ±0.50 


10205 


2 3 P 2 

Z 




10268.65 ±0.22 ±0.50 


10233 








10400 


3^ 


T(3S) 


10355.2 ±0.5 


10450 


2 1 P> 2 






10505 


2 3 D X 






10495 


2 3 D 2 






10502 


2 3 P> 3 






10511 


3^1 






10580 


3 3 P 






10562 


3 3 Pi 






10571 


3 3 P 2 






10589 


4% 






10716 


4 3 ^ 


T{4S) 


10579.4 ±1.2 


10764 
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